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Abstract 

With the help of the generalized characteristics(GC) of the first order partial differential equa- 
tions (PDE) we calculate the differential equation system of characteristics of the homogenous 
magneto hydro dynamical equations(MHD). 
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I. INTRODUCTION 



The investigation of the magneto hydrodynamical(MHD) equations is crucial importance 
to understand plasma instabilities in the future fusion facilities like ITER. Nowadays large 
efforts are made to numerically solve the complete MHD equations on a torus or on a more 
realistic stellarator geometry. l| However, simplified models of the complete MHD equation 
may help us to get a deeper understanding physical mechanism of plasma or let more insight 
into the structure of magnetized fluids. In the following study we shortly introduce the 
mathematical formalism of the generalized characteristics of first order partial differential 
equations and apply this theory to the full three dimensional MHD equations. 



II. THEORY 



It is well known from the theory of first order linear PDEs that along the characteristic 
line the original PDE becomes an ordinary differential equation(ODE). Therefore, we can 
make an important observation, along the characteristics the solution is constant. For 
linear equations the characteristic curve is a line. For one space and time dimension the 
characteristic lines give us a qualitative picture about the solution of the equation on the 
plain. The proof of this statement can be found in any textbook Q 

It the following we briefly introduce the mathematics of the generalized characteristics of 
the fist order PDEs. Let's consider the following first order PDE, which can be non-linear 
as well: 



F(x,t,u(x,t),p,q) = (1) 

we use the standard notation of p = du(x,t)/dx, q = du(x,t)/dt. According to the book of 
Melikyan Q] the differential equation system of the characteristics is the following: 

x = F p , i = F q , (2) 
u(x,t) = p ■ F p + q ■ F q , (3) 
p=-F x +p-F Ul q = -F t + q-F u (4) 

where . = 4-. To avoid further confusions and misunderstanding we use the more detailed 
notation: 
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dx _ dF dt dF 

9m dF du OF 

~& = di' d{§) + 3(f)' 
dp dF du dF dq dF du dF 

dr dx dx du ' dr dt dt du 

(8) 

These equations can be easily generalized to a higher dimensional equation system % = 
l...n as well. 

F i (x,u,z,t,Ui(x,u,z,t),pi,q i ,r i ,Si) = (9) 

with the notation of Pi = dui(x,y, z,t)/dx , % = dui(x,y, z,t)/dy, = dui(x,y,z,t)/dz, 
Si = dui(x,y, z,t)/dt. with the above given detailed notation 

dx dF,-, 



dr d{^) 
dy dFj 

dr d(^) 
dz dF 



dr d(f^) 
dt dF 



dr 

dui dui dFi dui dFi dui dFi dui dFi 



dr dx d(^) dy d(^) dz d(^) dt d(^) 

dpi dFi dui dFi 

dr dx dx du^ 

dqi dFi dui dFi 

dr dy dy dui 

dri dFi dui dFi 

dr dz dz dui 

dsi dFi dui dFi 



(10) 



dr dt dt dui 

These kind of equations can be applied to any kind of hyperbolic equation systems, like gas 
dynamical problems multi-phase flows [s| or to MHD. 
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The governing equations for an ideal, non-relativistic compressible plasma may be written 
in different forms if the following assumptions hold: 



A 



<1, 



<1, 



«1, 



pVL 



< 1 



(11) 
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where p,v,r and L are, respectively, characteristic density, speed, time and length scales for 
the problem, c is the speed of light, and e and a is the dielectric constant and conductivity of 
the fluid. In conservative variables, the governing equations, which is a combination of Euler 
equations of gas dynamics and the Maxwell equations of electromagnetics, is the following: 
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pvv + I (tt + ^p) 
vB - Bv 

V (E + tt + ^p) v - B(v • B) J 

where I is the 3x3 identity matrix, p is the density, v is the velocity, tt is the pressure (to 
avoid confusion with the notation used for partial differential equations p = du(x, t)/dx), B 
is the magnetic field, and E is the energy, defined as: 



E 



71 



+ P- 



v • V 



B B 



(13) 



7 - 1 ' r 2 2 

Solution of these equations can help to understand a number of problems governed by fluid- 
dynamics and electromagnetic effects. 

Given the following primitive variables w = (p, v x , v y , v z , B x , B y , B z , tt), the MHD equa- 
tions ffl2l may be written in quasi-linear form as: 



dw ^ w C^ W 

dt =dx = dy = dz 

where A, B_, C_ are 8x8 matrices. 
With a 
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where 7 is the compressibility of the fluid. Matrix A is singular - the fifth is zero, leading 
to a zero eigenvalue after the diagonalization. 



A 6 



A 7 



Vx + 



2p 



\ A - 7, 4- Bx \ A - 
, A 3 — V x H — , A 4 — 

V P 


B x 

v x — — 

VP 


B 2 + 7 7T+ v/(B 2 + 7vr) 2 


- 477TS 2 


B 2 + 7 7T+ v / (B 2 + 77t) 2 


- 477T5 2 


B 2 + 7 7T- v/(B 2 + 7tt) 2 


- 477T-B 2 


B 2 + 7 7t - v/(B 2 + 7 tt) 2 


4 7 7Tfi 2 



(16) 
(17) 
(18) 
(19) 



The eigenvalues of the system are well known, and they correspond to: 

• one entropy wave X A traveling with speed v x , 

• two Alfven waves A^ 3 traveling with speed v x ± c c where c a = ^= is the Alfven speed, 

• four magneto- acoustic waves ( \ A ...\ A ) 

For completeness we calculate and present the eigenvalues for the y and z direction too. 
Matrices B, C and the corresponding eigenvalues are very similar to A and have the same 
structure. 
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Xf=Vy + 



\ if 

A 6 = Vy ~ 
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The matrix form of the MHD equations (JT3J) reads as follows: 

I- 

dv x dv x dv x dv x 

-rrr + v x-^~ + v v^r + v z~k- = 
at ax ay az 



dt \ yfp J dx \ yfp J dy \ yfp ) dz 
dv z ( B x \dv z f B y \ dv z ( B z \ dv z 

+ ~w P \I B2 + 77r + V^ 2 + 77r)2 - 477r ^ ) ^ + 



9z 



as, 
at 



+ v x 



\=\jw + 7 7T + v/(B 2 + 7 tt) 2 - 4 7 ^f J 
-I=^/b 2 + 7 vr + v^B 2 " -f 7 vr) 2 - ^nB^j ^ + 

v„ - ^B 2 + 7 tt + y/ (B 2 + 7 tt) 2 - 4^5^ ^ + 

„ 



ae, 
at 



^\/B 2 + 7 7T + v/(B 2 + 7 7T) 2 - 4 7 7T J B 2 J 

-^=^/B 2 + 7 vr - v/(B 2 + 7 tt) 2 - 4 7 vrS 2 ^ ^ + 



Vx + 



I / \ dB 

^ + ^=\/B 2 + 7 7T- v /(B 2 + 7 7r) 2 - 4 7 7r J B 2 ^ = 



2p v ' v v ' y 1 Z J dz 

^r+(v x - -I^B 2 + 7 tt - v/(B 2 + 7 vr) 2 -~4^ef) ^ + 



at 



(30) 



' " Wp V B2 + 77r " V /(B2 + 77r)2 " 477r ^ 2 ) f + 

U - -^=^/B 2 + 7 vr - v/(B 2 + 7 vr) 2 - 4 7 ttB^ ^ = 

Now applying the equations of ({TO]) , after a tedious derivation we may get the complete 
equation system of the generalized characteristics. 

For the first equation of (|30|) §f = the equation system of characteristics is trivial 
Jr = 1, ^ = s where r is the parameter of the characteristic curve. It is clear that time can 
be used as a natural parameter too, so the equation of the first characteristics is — Si. 
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For the second equation of (130|) the system becomes much more complicated: 
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to avoid too lengthy formulas we used the standard eigenvalue notations 
(A^, j — A,B,C k — 5, 6, 7, 8) for the last four magneto-acoustic waves. The complete 
system of equations for the characteristic curves are 8 x 8 = 64 equations, unfortunately 
the author has problems to interpret the results in the present form. We can restrict the 
motion to reduce the complexity of the dynamics. The simplest case is a one dimensional 
motion in a one dimensional magnetic field which means a system with 4 PDEs. A possible 
stability analysis of a restricted motion may give new insight into the dynamics of the 
MHD equations. We can imagine that even more complex analysis like the Panleve test or 
some kind of stability analysis can be applied to the presented differential equation system, 
unfortunately such investigation are out of the skill of the author. We can only hope that our 
work may give some impetus and motivate some reader to start such studies and investigate 
the equations presented above. 

III. SUMMARY 

We briefly presented how the generalized characteristics of the homogenous magneto 
hydrodynamical equations can be derived. Unfortunately, no deeper e.g. geometrical in- 
terpretation of the equation could be found, other detailed analysis is also lacking or may 
remain to the reader. 
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